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Using low-temperature scanning tunneling spectroscopy applied to the Cs-induced two- 
dimensional electron system (2DES) on p-type InSb(llO), we probe electron-electron interaction 
effects in the quantum Hall regime. The 2DES is decoupled from bulk states and exhibits spreading 
resistance within the insulating quantum Hall phases. In quantitative agreement with calculations 
we find an exchange enhancement of the spin splitting. Moreover, we observe that both the spatially 
averaged as well as the local density of states feature a characteristic Coulomb gap at the Fermi 
level. These results show that electron-electron interaction can be probed down to a resolution 
below all relevant length scales. 

PACS numbers: 73.21. Fg, 73.61.Ey, 71.70.Gm, 68.37.Ef 



Quantum Hall (QH) physics [1] is a paradigm for the 
study of interacting quantum systems [2]. In this respect, 
the III-V semiconductors are the very mature materials, 
although graphene catches up [3]. The most intriguing 
QH phases are driven by electron-electron (e-e) interac- 
tion [4, 5], which, however, is screened by nearby gates 
and competes with disorder. Thus, a central challenge 
towards a microscopic investigation of QH physics domi- 
nated by e-e interaction is to provide a sufficiently clean 
and electrically decoupled system probed down to the 
relevant length scales, most notably the magnetic length 
Ib = v h/(eB) ~ 10 nm (6 T). Scanning tunneling spec- 
troscopy (STS) achieves the required nm resolution. Ap- 
plying STS to an adsorbate induced 2D electron system 
(2DES) [6, 7], some of us have shown that the states re- 
sponsible for the integer QH transitions can indeed be 
probed with nm resolution [8, 9]. Theoretical analysis 
[10] of these data and similar experimental results for 
graphite and graphene [11] have also been published. 

Here, we modify the 2DES in order to fully decouple 
it from the substrate and to reduce the disorder. This 
allows to probe e-e interaction effects. In particular, we 
observe an exchange enhancement (EE) of the spin split- 
ting at odd filling factors in quantitative agreement with 
a parameter-free calculation. Moreover, we measure a 
Coulomb gap in the spatially averaged density of states 
(DOS) at the Fermi level Ep. This Coulomb suppression 
is in quantitative agreement with predictions for localized 
systems [12-14]. Interestingly, we find a similar suppres- 
sion in the local DOS (LDOS), which is probably caused 
by fluctuation effects. Observing these hallmarks of the 
e-e interaction in STS is a crucial step towards a direct 
imaging of intriguing QH states such as stripe, bubble or 
fractional QH [4, 15] phases. 

The home-built scanning tunneling microscope oper- 
ates at T = 5 K in ultra-high vacuum (UHV) [16]. The 
dl/dV curves representing the LDOS of the sample are 



measured by lock-in technique at constant tip-surface 
distance stabilized at current 7 s t a b and voltage T4tab- A 
modulation voltage V mo d is used to detect dl/dV while 
ramping the sample voltage V. 

The 2DES was prepared in UHV by cleavage of a p- 
type InSb single crystal (Na = 1.1 X 10 21 m~ 3 ) and sub- 
sequent Cs adsorption of 1.1 % of a monolayer (3.7 x 
10 16 m~ 2 ) onto the cooled (110) surface [17]. The single 
Cs atoms act as surface donors, which bend the bands 
downwards and induce a 2DES [6, 9, 18]. Fig. 1(a) shows 
the corresponding band bending in the near-surface re- 
gion as calculated by a Poisson solver. This leads to a 
2DES with density iV s ~ 2.7 x 10 16 m" 2 . Note that 
the band bending reaches deep into the bulk leading 
to a decoupling of the confined states of the 2DES at 
Ei, i £ No, from the partly empty bulk valence band 
(BVB) 600 nm apart. Indeed, the spatially averaged 
dl/dV curve (I) of Fig. 1(b), measured without contact- 
ing the 2DES directly, does not exhibit any signature of 
the 2DES, but only an increase in dl/dV close to the on- 
set of the bulk conduction band (+200 mV) and the sur- 
face valence band (—400 mV). The tunneling path from 
the 2DES to the BVB is blocked. This is in contrast to 
measurements using n-type [7, 9, 19] and p-type samples 
with higher doping [20, 21], always exhibiting a step like 
increase in spatially averaged dl/dV curves close to the 
calculated Ei. If our 2DES is additionally contacted by 
an Ag stripe running perpendicular to the cleavage plane 
[22], it exhibits two steps close to the calculated Eq and 
Ei as visible in curve (II) of Fig. 1(b). 

Applying a magnetic field B perpendicular to the 
2DES results in peaks corresponding to Landau levels 
(LLs) and spin levels of the 2DES [curve (III), Fig. 1(b)]. 
Their distance is in accordance with the effective mass 
m* and ^-factor g* of the InSb conduction band [17]. 
The width of the peaks at lower energy is AE ~ 16 meV 
which is caused by the potential disorder, mainly given 
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FIG. 2. (color online), (a) dl /dV spectra (grayscale) recorded 
at the same lateral position at different B fields as indicated, 
/stab is increased for each B from 100 pA to 2000 pA (left to 
right). Vstah = 300 mV, V mo d = 1.6 mVrms. (b) Measured 
lowest LL positions (spin up t; spin down 4-) at B = 7 T 
(symbols) in comparison with calculated LL positions (line). 
Note that the current I at the peak position and not /stab is 
used as x axis. 



FIG. 1. (color online), (a) Calculated band bending us- 
ing ID Poisson equation; first two 2D subband energies 
Eo = —118 meV, E\ — —65 meV as calculated using the tri- 
angular well approximation [17, 23] are marked; correspond- 
ing electron distribution is shown in [17]. (b) Spatially aver- 
aged dl/dV curves across an area Aa v '- (I): without contacted 
2DES, (II), (III): with contacted 2DES at B = T, 7 T as 
marked; gray area: bulk band gap of p-InSb [24]; Eo, E\: sub- 
band energies from experiment (solid) and calculation (broken 
lines); arrows in (III) mark spin split LLs. (I) VJtab = 400 mV, 
/stab = 100 pA, V mod = 1 mV rms , Aa v = 200 x 160 ran 2 ; 
(II) Ktab = 300 mV, l s tab = 50 pA, V mod = 3 mV rms , 
A Av = 100 x 100 nm 2 ; (III) T4tab = 300 mV, / st ab = 200 pA, 
Vmod = 0.4 mVrms, A Av = 300 x 300 nm 2 . 



by the dopants of the substrate. The Cs atoms, which 
are ionized by only 30 % [25] , have a minor effect [7] . 

Further evidence for the electrical decoupling of the 
2DES from the BVB is presented in Fig. 2, where dl/dV 
curves at increasing I s tab are shown for different B. All 
curves are measured at the same lateral position. With 
increasing B, pairs of lines corresponding to spin split 
LLs appear. For B = T, increasing I a tab does not 
change the dl/dV spectra. At higher B the spectra are 
spread in V with increasing / s tab- At 7 T, the spin split- 
ting of the lowest LL is increased by 13 % and the LL dis- 
tance is increased by 18 %. The spreading is symmetric 
around V = mV, i.e. around E F . It is attributed to the 
increased localization of electrons with growing B lead- 
ing to a decrease of 2D conductivity [2]. The spreading 
cannot be explained by tip induced band bending (Stark 
effect) [26], which might increase with B due to a re- 
duced screening of the 2DES. Poisson calculations reveal 
that the spreading at largest tip-surface distance must be 
much larger than the spreading induced by the change in 
tip-surface distance in contrast to experiment. Instead, 
the spreading is quantitatively reproduced by assuming 
a thermally activated nearest neighbor hopping of the lo- 
calized electrons within the 2DES from or towards the 
tip. The model uses barrier heights and next-neighbor 



distances of valleys as determined from spatially resolved 
dl/dV data [17] and assumes a reasonable attempt fre- 
quency of vq = 10 13 Hz [27]. The resulting peak positions 
in comparison with the experimental data for B = 7 T 
are shown in Fig. 2(b). The excellent agreement strongly 
supports our assumption that the current indeed flows 
along the 2DES exhibiting reduced conductivity with in- 
creasing B. 

The surface 2DES, thus, is occupied, exhibits Lan- 
dau as well as spin quantization, has moderate disorder, 
and is decoupled from the bulk electrons of InSb. More- 
over, the center of mass of the 2DES is 8 nm below the 
surface and, thus, sufficiently far from the metallic tip 
to prevent complete screening. These are the require- 
ments to observe e-e interaction effects within the QH 
regime. One such effect is the EE of the spin split- 
ting. Loosely speaking the effective repulsion between 
electrons with parallel spins is smaller than the one for 
antiparallel spins. This eventually leads to an increase 
of the spin splitting energy Ess a t odd filling factors 
[28]. Fig. 3(a) shows dl/dV spectra taken at a fixed 
position while ramping the magnetic field B providing 
the so-called Landau fan. Less than 10 % of the fanning 
is caused by the spreading resistance described above. 
Varying B, the conductance lines are wavy and deviate 
from _E" ± = Ei + fuj c (n + \) ± \g* ^B with subband 
index i, LL index n £ No, spin index ±, cyclotron fre- 
quency uj c = eB/m*, and Bohr magneton fi^- One ob- 
vious reason for waviness is a shift of E-p with magnetic 
field taking place once the increasing degeneracy of a LL 
favors a transition to the next LL. More importantly, g* 
is filling factor dependent due to the EE. To analyze this 
in more detail, we concentrate on the lowest LL around 
— 120 mV, which gives the highest accuracy in determin- 
ing Ess- We adapted two Gaussians for all 386 spectra 
between 3.5 and 6.1 T [29]. The fits are good as can be 
seen in the inset of Fig. 3(b) and by the confidence value 
of R 2 = 0.94 (0.97 above 5 T). The error for the resulting 
Ess is about 0.2 meV. Ess(B) is shown in Fig. 3(b) in 
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FIG. 3. (color online), (a) Landau fan showing dl/dV 
as grayscale; B field ramped downwards, V s tab = 300 mV, 
4tab = 400 pA, Vmod = 1.6 mV rms . (b) Spin splitting en- 
ergy Ess of lowest LL extracted by Gaussian fits as shown 
in the inset. Straight line marks .Egg = |<?*|mb^- (c) Devi- 
ation A from the linear fit in (b). The inner bright line is 
smoothed. Vertical bars show the calculated values of EE be- 
tween neighboring odd and even filling factors. Local filling 
factors v oc 1/B are matched by counting spin levels below 
E-p (compare Fig. 4(b)) and verifying their B dependence. 



comparison to a straight line corresponding to ordinary 
Zeeman splitting of \g*\n&B with \g*\ = 42. Fig. 3(c) 
shows the deviation A(B) from the straight line. It os- 
cillates around meV with maxima (minima) around odd 
(even) filling factors as expected for EE [30] . The not ex- 
pected negative values of A(B) are probably caused by 
slight deviations from a spin splitting linear in B due to 
either increasing spreading with B, which leads to super- 
linearity, or nonparabolicity of InSb leading to a smooth 
decrease of g*, thus, supralincarity. However, both effects 
cannot explain oscillations in A(B). One could imagine 
that spreading depends also on filling factor being largest 
at even ones, but that would lead to an oscillation with 
maxima at even filling factor. 

Moreover, the amplitude of the A(B) oscillation is 
about 0.7 meV in excellent agreement with theoretical 
estimates for EE [vertical bars in Fig. 3(c)]. They are 
obtained by treating the Coulomb interaction using a 
random phase approximation. This is well justified since 
the subband electron density TVo is large compared to the 
scale set by the Bohr radius [17, 28, 31]. We performed 
the calculation using m* = 0.02mo and g* = —42 but em- 
phasize that the results barely change if these or other 
system parameters are varied within reasonable limits 



\ B = 7T 









-20 





V(mV) 



20 




CO 

iT 0.5 



peak 
max. 



B = 5.56 T, v = 9.6 

B = 5.50 T, v = 9.7 T 



-10 

V(mV) 



10 




-100 -50 
l/(mV) 



FIG. 4. (color online), (a) Zoom of Fig. l(b)(III) (thick line), 
additional Gaussian fits for the LL at Ef (thin lines), and a 
plot of the difference (medium, blue line) in comparison with 
the expected bare Coulomb gap of eq. 1 (V-shaped line) and 
the Coulomb gap taking finite temperature, screening and 
energy resolution into account (dotted line); absolute scale of 
DOS is deduced by fitting the lowest spin-split LL to a two- 
peak Gaussian and matching the integral of one Gaussian to 
the known level degeneracy eB/(2nh). (b) dl/dV spectra 
taken at the same position at B as marked. Arrows of same 
color follow the same peak or minimum across Ef', a movie 
of the data is available in [17] (c) dl jdV value of the marked 
peak and minimum in (b) as a function of energy (voltage). 



(e.g. less than 1 % for g* = —38). Thus, magnitude 
and oscillation phase of A(B) compare favorably with a 
parameter-free calculation of EE. This implies that the 
short-ranged e-e interaction effect EE can be probed by 
STS. 

For localized electrons interacting via the long-ranged 
part of the Coulomb repulsion, the averaged tunneling 
DOS is expected to show a gap at Ef [12-14]. For a 
2DES with unscreened repulsion at T = K, a qualita- 
tive analysis gives [13, 14] 



D (E)^-^1\E-Ef\ 



(1) 



More elaborate analytical and numerical results leave 
no doubt about the existence of a Coulomb gap while 
the exact shape remains controversial [14, 32]. This is 
due to the underlying (spin-)glass physics [32] known 
to be notoriously complex. A linear Coulomb gap was 
deduced from various experiments [33]. In our case, 
where the ratio between disorder and e-e interaction is 
R = AE/[(e 2 y/N~ s )/(4irs r s a )} ~ 1.1, we also find a dip 
in the DOS at Ef- Fig. 4(a) shows the spatially aver- 
aged dl/dV curve (thick line) at B = 7 T. Instead of 
a peak at Ef, one observes a double-peak with a mini- 
mum at mV. The sum of two identical Gaussian peaks 
(thin lines) — mimicking the two spin levels of this par- 
ticular LL (sec [17]) — matches the measured DOS except 
of a suppression at Ef- Taking the difference between 
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measured DOS and the sum of the two Gaussians elim- 
inates all single-particle effects leaving only the dip at 
Ep (medium line). If we modify Eq. (1) to account for 
finite temperature and screening effects [14] as well as 
for the energy resolution of our experiment of 1.6 meV 
[17], we obtain the dotted curve in Fig. 4(a). It shows 
excellent agreement with the measured dip. The screen- 
ing is taken to be caused by the STM tip being 8.6 nm 
away from the center of mass of the 2DES [17]. Note 
that we observe the gap even around the critical state, 
i.e., close to half filling of a spin-polarized LL, which is 
consistent with numerical studies [34]. The facts that 
we do not observe the dip at Ep without localization (at 
B = 0T) and that we can reproduce it by a reasonable, 
parameter-free calculation strongly suggests that we ob- 
serve the Coulomb gap. We can rule out inelastic exci- 
tations as a cause which would lead to much larger half 
widths of the gap (optical phonons: 22 meV, plasmons: 
60 meV, spin excitations: 18 meV) and we are not aware 
of any many-particle mechanism besides the long-ranged 
Coulomb repulsion of localized electrons leading to a gap 
with the observed characteristics. 

Surprisingly, a Coulomb gap — although typically 
thought of being a phenomenon related to disorder aver- 
aging or spatial averaging — is also observed in the local 
DOS [35]. The intensity of a particular LDOS peak is 
suppressed when moved through E-p by increasing B. 
Fig. 4(b) shows corresponding dl/dV curves at fixed 
position. The upper arrows follow a single spin level 
as it crosses Ep and the peak intensity is plotted in 
Fig. 4(c). A minimum intensity is observed exactly at 
Ep (B = 5.62 T) where the peak is suppressed by 48 % 
(suppression in averaged DOS: 33 %). The same kind 
of dl/dV suppression is found for the minimum between 
LLs, which is marked by the lower lying arrows in (b) and 
plotted in (c), too. A dl/dV suppression at Ep is also 
found for fixed B, if different positions are probed within 
the potential landscape [17]. The finding of Coulomb 
suppression in the LDOS requires further studies and 
might be related to Coulomb glass dynamics [36]. 

In summary, we have shown that low-temperature STS 
is able to detect e-e interaction in QH samples down to a 
resolution below all relevant length scales. We have found 
an exchange enhancement (EE) of the spin splitting at 
odd fillings and a Coulomb suppression of the averaged 
as well as of the local DOS at Ep. The EE is in quanti- 
tative agreement with a well justified theory, while, due 
to the less clear status of theory, the comparison for the 
Coulomb gap is with calculations based on qualitative 
arguments only. No well-developed theory for the LDOS 
exists and we conjecture that the Coulomb gap in LDOS 
is related to (spin-) glass physics. 

We acknowledge support by the DFG (MO 858/11-1). 




FIG. SI. STM image of InSb(llO) with 1.1 % Cs coverage; 

V = 300 mV, I = 30 pA, 300 nm x 300 nm. The inset shows 
an area of 50 nm x 50 nm from a separate measurement with 

V = 300 mV, I = 50 pA. 



SUPPLEMENTARY INFORMATION 
EXPERIMENTAL DETAILS 

For sample preparation, we glued a p-InSb single crys- 
tal of size 3 mm x 3 mm x 3 cm to a sample holder using 
silver epoxy. A silver epoxy line drawn at one side of the 
crystal from the sample holder towards the surface serves 
as a direct electrical contact to the 2DES. At the opposite 
side of the epoxy line, the crystal was cut for controlled 
cleavage about 1 mm in depth. This cut is parallel to the 
sample holder surface and located about 1-2 mm above 
the sample holder. After bake-out within the load-lock 
of our ultrahigh vacuum system, the crystal was cleaved 
at room temperature and a pressure of p ~ 10~ 8 Pa 
exhibiting the (110) surface. Then, we transferred the 
crystal into the pre-cooled scanning tunneling microscope 
(STM). For this purpose, the STM was lifted from the 
bath cryostat to the transfer chamber without breaking 
the vacuum [16]. Cesium from a well-outgasscd dispenser 
(SAES Getters) was evaporated onto the cold crystal sur- 
face. During and after the adsorption process the sample 
was held below 50 K to prevent diffusion induced Cs clus- 
tering [21]. Fig. SI shows a resulting STM image. The 
image area is atomically flat and each bright dot corre- 
sponds to a single cesium atom or, partly, to a Cs dimcr 
which both act as donors. 3300 adsor bates are visible 
which corresponds to a coverage of 1.1 % per InSb(llO) 
surface unit cell. 
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FIG. S2. Spatially averaged dl /dV curve (black line) at B = 
7 T, V stab = 300 mV, / stab = 200 pA, V mod = 0.4 mV rms , 
Aav = 300 nm x 300 nm, 20 x 20 spectra; yellow areas show 
the individual Gaussian fit curves of the spin split Landau 
levels marked by arrows each; red line marks the resulting fit 
of the dl/dV curve gained by adding up the Gaussians; peak 
energies of the Gaussians are labeled above each peak in meV. 

The STM tip was etched from a tungsten wire outside 
of the vacuum system and prepared within the STM by 
field emission and consecutive voltage pulses on a W(110) 
crystal as well as on the InSb(llO) surface itself. After 
such voltage pulses on InSb, the scan area had to be 
changed using the translation stage of the STM. We used 
tungsten because it enabled us to give good results on a 
similar sample system not exhibiting tip-induced band 
bending [21]. 



LANDAU LEVELS AND POTENTIAL DISORDER 

At B = 7 T, we observe the spin split Landau levels 
of the first subband as shown in Fig. 1(b) of the main 
article. The same measurement is shown in Fig. S2 with 
the corresponding Gaussian fits. From the positions of 
the fits we extract a spin splitting for the first Landau 
level of 16.7 meV = |<j*|/zb-B giving an effective Lande 
g-factor of \g*\ = 41. From the mean difference of the 
first two Landau levels we get an effective mass of m* = 
(heB)/(37 meV) = 0.022 x mo (mo : free electron mass). 
These values are in good agreement with previous STS 
measurements on InSb [9, 21] and differ slightly from the 
known low temperature values for InSb at the conduction 
band minimum m^/mo = 0.0135 [24] and g = —51 [37] 
due to the well-known nonparabolicity of the conduction 
band of InSb. 

The full width at half maximum (FWHM) of the fits 
decreases with Landau level index n, being 16.0 mV for 
n = 0, 12.7 mV for n = 1 and 12.1 mV for the first spin 
level of n = 2. This effect is due to the decreased sensitiv- 
ity of the Landau level drift states to the potential land- 
scape, which can only probe potential fluctuations down 
to length scales of r c = y/(2n+ l)h/(\e\B) [9]. Single 
spot dl/dV spectra only show a width of about 11 mV 
for n = 0. From this, we can estimate the width of the po- 




FIG. S3. dI/dV(x,y) image at V = -132 mV and B = 
7 T, V stab = 300 mV, 7 stab = 200 pA, V mod = 1.6 mV rms , 
300 nm x 300 nm. High dl /dV intensity marks potential pits. 



tential disorder to be Adis — Vl6 2 — ll 2 meV = 12 meV. 

The main cause for the potential disorder are the 
charged acceptors within the subband layer. The po- 
tential disorder of a highly doped p-InSb(llO) sample 
covered with Cs has been larger with Adis — 20 meV 
[21]. A previously analyzed 2DES prepared by Cs cov- 
erage on n-InSb(llO) showed a potential disorder Adi S 
comparable to our system with Adis — 10-15 meV, de- 
duced from the lowest spin level at T = 300 mK and 
B = 12 T [9]. The disorder of a 2DES created by Fe 
adsorption on n-InAs(llO) showed a potential disorder 
of about Adis — 20 meV [7, 8], which was determined by 
two independent methods: an STM image of the poten- 
tial landscape was generated by probing the surface with 
a tip-induced quantum dot state [7] and FWHM of Lan- 
dau level peaks measured at T = 6 K and B = 6 T were 
evaluated [8]. However, these measurements are not spin 
resolved due to the lower g-factor of In As, so a direct 
comparison of Adi s is difficult. 

The lateral distribution of the potential disorder can 
be imaged by scanning tunneling spectroscopy (STS) in 
real space as is shown in Fig. S3. The voltage of this 
dl/dV image corresponds to the onset of the first Lan- 
dau level and, thus, a high dl/dV intensity exhibits the 
potential pits. There are roughly 40 local potential min- 
ima within the imaged area, giving an average distance 
of ^(300 nm) 2 /40 ~ 50 nm in accordance with visual 
inspection. 
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BAND BENDING AND SUBBAND 
CALCULATION 

It is well-known that submonolaycr amounts of atoms 
adsorbed on an InSb surface lead to a large band bending 
near the surface. The adsorbed atoms act as donors to 
the underlying substrate with a donor level in the con- 
duction band. For cesium on n-InSb(HO), this donor 
level has been determined by photoelectron spectroscopy 
to be 290 meV above the conduction band minimum [18]. 
In our case of p-doped InSb, where the bulk Fermi level 
Ep lies at the valence band edge, we have to add the 
low temperature band gap of E gap = 235 meV [24] to 
get an estimate of the total band shift at the surface 
of Vbb = 525 meV. The corresponding band bending 
leads to electrons within an inversion layer which form 
two-dimensional subbands. The curvature of the poten- 
tial inside the semiconductor, therefore, depends on the 
charge distribution within the inversion layer and on the 
density of charged acceptors. 

The sample we used has an acceptor concentration of 
TVa = 1.1 x 10 21 m~ 3 , which is orders of magnitude lower 
than the inversion layer electron concentration. To get an 
estimate of the expected subband energies without hav- 
ing to solve the Poisson equation and the full Hamiltonian 
including the nonparabolicity of InSb self-consistently, we 
approximate the inversion layer distribution by the 3D 
bulk density of states: 

Dm(E) = &%MiVE. (SI) 

The energy dependence of the effective mass m* due 
to the nonparabolicity is approximated by m*(E) = 
m^il + 2E/E g3iP ) [38] with the known effective mass at 
the conduction band minimum raj/mo = 0.0135 [24]. 
This results in a band bending as shown in Fig. 1(a) of 
the main article with an inversion layer electron concen- 
tration of N s = 2.7 x 10 16 m~ 2 . 

To obtain the subband positions Ei, i G No within 
this potential, we use the triangular well approximation 
[39] with electric fields approximated by the mean slope 
of the potential up to Ei and an effective mass using the 
approximation m*(E) = m* (l + 2{E,/3 + E ll )/E gllp ) [38]. 
Here, the energy E splits up into the in-plane energy E\\ 
and the energy Ei at the onset of the subbands. The 
subband positions are marked in Fig. 1(a) by horizontal 
lines. The density of subbands above the Fermi level 
is relatively large because of the flat band bending at 
a depth z beyond the inversion layer, where only the 
acceptor density can screen the electric field. The first 
two subbands E and Ex are in good agreement with the 
experiment as shown in Fig. 1(b). 

The 2D electron concentrations of the first two sub- 
bands are N Q = 1-3 x 10 16 m~ 2 and N\ = 0.7 x 10 16 m" 2 . 
Their electron distributions along z arc plotted in Fig. S4. 
From these distributions, the mean depth of the subband 
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FIG. S4. Zoom of Fig. 1(a): Calculated band bending (black 
line) and the first two 2D subbands with their electron distri- 
bution curves |i/>i(,z)| 2 (yellow areas). The mean depths (z) of 
each distribution is indicated by the dotted lines and marked 
in nm. 

electrons below the surface is deduced to be 7.9 nm for 
the first subband and 15.2 nm for the second subband. 

The averaged dl/dV curve (curve (II) of Fig. 1(b)) 
exhibits rounded steps close to the energies calculated 
by the Poisson-Schrodinger equation. The width of these 
steps is caused by disorder. Indeed, spatially resolved 
dl/dV images as shown e.g. in Fig. 3 of Ref. [7] reveal 
quantum dot like states at low energy which percolate 
close to the onset of the plateau. The decreasing height 
of the step with increasing subband index is consistently 
observed for different 2DES [7, 20, 21] and is attributed 
to the fact that higher subbands are reaching deeper into 
the crystal and, thus, have a lower LDOS at the surface. 
This is clearly visible in Fig. S4. 



COULOMB GAP AND ENERGY RESOLUTION 

In order to account for the finite temperature T = 5 K 
and for screening effects by the metallic STM tip, we 
have modified the V-shaped DOS of Eq. (1) of the main 
article called Dq(E) by a finite DOS at Ef according to 
Rcf. 14: 

D mod (E F ) « 0.085^ + . 86 (^) 2 fc B T 
e d e 

with d being the distance of a metallic gate to the cen- 
ter of the 2DES and fee being the Boltzmann constant. 
In our system, d can be estimated by the mean depth 
of the first subband as calculated above and the dis- 
tance of the STM tip to the surface. The distance of 
the tip to the surface is estimated by fitting I(z) mea- 
surements to an exponential tunneling decay and extrap- 
olating the z = offset to the quantum of conductance 
Go = 2e 2 /h with electron charge e and Planck's con- 
stant h [40]. From this, we get a tip distance of about 
700 pm at the stabilizing current of this measurement 
/stab = 400 pA. This results in d = 8.6 nm. With 
e r = 16.8 [41], we get D mod (E F ) = 1.66 x 10 17 cV^nT 2 . 
For the sake of simplicity, we simply cut off the V- 



shaped Coulomb gap Dq(E) by this value leading to 
D(E) = max{D (E),D mod (E F )}. 

The Coulomb gap D(E) is finally broadened by the 
energy resolution of the STS experiment being 8E ~ 
^(3fc B T) 2 + (2.5Knod) 2 = 4.2 meV [42, 43] with V mod = 
1.6 mV rms and T = 5 K. In order to calculate the re- 
sulting dl/dV curve, we approximate a Gaussian energy 
broadening with SE being twice its standard deviation. 
This Gaussian is folded with the calculated D(E) result- 
ing in the red dotted curve of Fig. 4(a) of the main 
article. Notice that no fit parameter is involved within 
this calculation which excellently reproduces the mea- 
sured Coulomb gap. 



LOCAL COULOMB GAP WITHIN THE 
POTENTIAL LANDSCAPE 

The Coulomb gap in the DOS has been measured and 
discussed within the main article via spatially averaged 
dl/dV spectra. The main article and the attached video 
file shows, in addition, that single dl /dV peaks measured 
at a fixed position also exhibit a Coulomb gap, i.e. a re- 
duction of intensity while moved across Ep by B field. 
Such a Coulomb gap can also be identified at fixed mag- 
netic field in local spectra taken at different positions 
within the potential landscape. 

Fig. S5 shows a low energy dl/dV image marking the 
potential disorder of the same area as imaged in Fig. S3. 
Fig. S5 has lower spatial resolution, but higher energy 
resolution. The averged spectrum of all spectra gives the 
DOS of Fig. 1(b) (III) of the main article. Each pixel 
represents a spatial average of 12 neighboring spectra, 
thus exhibiting the averaged potential of an area of about 
50 nm 2 . The spectra corresponding to the marked posi- 
tions are shown in Fig. S6. 

Moving down in the potential landscape as indicated 
by the arrow, the peaks belonging to the three lowest 
Landau levels shift to lower energies as expected. Near 
the Fermi level, additional subbands lead to a more ir- 
regular pattern, but the Landau levels of the lowest sub- 
band still dominate the pattern. If observed carefully, 
all features shift to lower energies, except the suppres- 
sion caused by the Coulomb gap, which is pinned to 
V = mV. This suppression is visible as a dip at medium 
potential, where a spin state is exactly at E-p, and as 
kinks in the tails of the peaks at low and high poten- 
tial. This plot illustrates again the local observation of 
a Coulomb gap feature, now by averaging over distances 
belonging to a single potential minimum or maximum. 
Thus, the Coulomb gap appears to be ubiquitous even 
in local spectra, which is not expected from the typical 
qualitative arguments. 
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FIG. S5. dI/dV(x,y) image at V = -132 mV and B = 
7 T; the data belong to the spatially averaged dl/dV curve 
in Fig. 1(b) (III) of the main article; V sta b = 300 mV, 7 s tab = 
200 pA, Vmod = 0.4 mVrms, 300 nm x 300 nm, 20 x 20 spectra. 
High dl/dV intensity marks potential pits as in Fig. S3. The 
crosses and the arrow mark the positions and the sequence 
of the spectra shown in Fig. S6. 12 next neighboring spectra 
have been averaged for each spectrum in order to reduce noise. 
This leads to a spatial resolution of 60 nm as indicated by the 
circle. 




-150 -100 -50 50 

l/(mV) 

FIG. S6. dl/dV spectra recorded at different positions as 
marked in Fig. S5. The spectra are offset vertically. The 
Coulomb gap is visible as a dip or kink pinned at V = mV 
(E-p). Note that all other features except the Coulomb gap 
feature shift continuously due to the varying potential. 
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THEORETICAL ESTIMATE OF THE 
EXCHANGE ENHANCEMENT 



In this Section, we briefly sketch how to obtain a quan- 
titative estimate of the filling factor dependence of the ef- 
fective g-factor resulting from the Coulomb interaction. 
We start out by introducing the Hamiltonian describing 
a two-dimensional free electron gas subject to a perpen- 
dicular magnetic field B (pointing in z-dircction): 



Hn 



(S3) 



with Cnka being an annihilation operator in second quan- 
tization. The single-particle energies read 



flUJ r 



a eh B 



eB 



2 2mo c 



Oj. 



(S4) 



where the effective mass m* = 0.02mo and effective g- 
factor g* — —42 can both be extracted from our ex- 
perimental data. We have chosen the vector potential 
A{x) = Bxe y in the Landau gauge, and the n-th Landau 
level of electrons with spin a = ± is thus characterized 
by an additional quantum number k associated with the 
momentum operator p y of the y-direction. In contrast, 
the interaction is naturally expressed in terms of two- 



dimensional momentum quantum numbers k = (k x , k y ): 



U 



• 4 - c r H ■ ( S5 ) 

qtJx k 2 ~qo 2 k 2 a 2 k 1 a 1 



k 1 E 2 < 7 1 C7 2 

9/0 



Using cgs units, the Fourier transform of the Coulomb 
potential yields 



V(q) = F(q) 



2ne 2 
s T q 



(S6) 



with the form factor F(q) accounting for the finite thick- 
ness of the inversion layer [28] : 



F(q) 



3 

8x~ 



3 

8^2" 



1 

Ax 3 



3 / 487T7U ::: ( :-' 

e T h 2 



. 32 



The dielectric constant e r of InSb is given by e r 



N d ) 

(S7) 
16.8, 



and 7V S = 2.7 x 10 
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as well as iV d = 8 x 10 1 



are the total electron density in the inversion layer as well 
as the density of the ionized acceptors, respectively. 

In order to compute the self-energy £ of the many- 
particle problem posed by the Hamiltonian H = Hq + U , 
we employ the random phase approximation (RPA) [28] , 
which is reasonable for the problem at hand (at least for 
integer filling factors where the fractional quantum Hall 
effect is absent) since the density of the electrons occu- 
pying the lowest subband N « 25/a 2 , is large compared 
to the scale set by the Bohr radius ao = ti 2 £ 1 -/m*e 2 . The 
first (Fock exchange [30]) term of the RPA series reads 



S nifci,n 2 b( iw ) = ~ T zJ E ( n ikicr,n 3 k 3 a\U\n 3 k 3 o;n2k2a) (n 3 k 3 o-\Ga(iuj + iVL)\n 3 k 3 a) 

r»3&3 i& 



= -T5 mn2 S(h-k 2 ) J f^£ P -«3 (<z 2 / 2 /2)E 



V(q) 



iuj + itt- e nsa + fi 



(S8) 



where Qv(iu) is the nonintcracting Matsubara Green 
function related to Hq, I — \JhcjeB the magnetic length, 
and T = 5 K as well as fj, = vftw c j2 denote the temper- 
ature and the chemical potential, respectively. We have 
defined the following combination of associated Laguerre 
polynomials L™J _ni (x): 



Pn in2 (x) = (-1)" 1+ " 2 . 



'■(x)L^(x). (S9) 



In order to derive the second equality of Eq. (S8), 
one needs to evaluate the two-particle matrix element 
(nifcicr, n 3 k 3 a\U\n 3 k 3 a, n^kio) of the Coulomb interac- 
tion. This can be achieved by resorting to the momentum 



r 



representation of the Landau states \nka), 



(k x k y a\nka) = S(k — k y )(—i) 



i 22) (S10) 

■ e-* a / 3 ir„(ife), 



i!2" 



inserting unit operators, and eventually carrying out an 
integral over products of Hermite polynomials H n (x) 
[44]. 

The second term of the RPA series involves the bubble 
diagram depicted in Fig. S7. The corresponding analytic 
expression reads 
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TV(q)V(cf) 53Z)<*x^H^o(iA)|*xA,,tr) ((k x + q x ){k y + q y )a\g (iX + m)\(k' x + q' x ){k' y + q y )a) 



A- A/ 



with /(e) = l/{l + exp[(e-/i)/T]} being the Fermi func- 
tion. The second equality can again be established by 
inserting unit operators in order to calculate the (two- 
dimensional) momentum matrix elements of the nonin- 
teracting Green function. Since Eq. (Sll) is of the same 
(momentum-conserving) structure as the Coulomb po- 
tential, the geometric RPA series can be summed up in 
complete analogy with the well-known case where the 
magnetic field is absent [45]. The resulting self-energy is 
eventually given by Eq. (S8) with the bare V(q) replaced 
by 

v RPA (q, in) = , (si2) 

1 - V{q)Tl {q,iil) 

which wc moreover set to its zero-frequency value iQ = 
[28]. The remaining Matsubara sum in Eq. (S8) can then 
be carried out analytically, and the self-energy (which 
is frequency-independent, ruling out difficulties with the 
analytic continuation to the real axis from the beginning) 
can be computed with minor numerical effort. Surpris- 
ingly, it turns out that such a static screening approxi- 
mation is sufficient to reproduce the results of more elab- 
orate approaches [46] on a quantitative level, providing 
the a posteriori motivation to solely stick to this scheme. 

Since the self-energy is diagonal [47] w.r.t. n and k, 
the magnitude of the oscillations shown in Fig. 4(c) of 
the main text is given by 

^) = (st =0 -Sl =0 )^-(st =0 -Sl =0 )^_ i , (S13) 

where the two bracketed terms are to be computed at 
filling factors v and v — 1, respectively. The result is 
(5(11) = 0.76 meV, <5(15) = 0.55 meV in nice agreement 
with the experimental data. We emphasize that these 
values only change slightly when the physical system pa- 
rameters (such as the dielectric constant, the tempera- 
ture, the effective (/-factor, or the effective mass) are var- 
ied within reasonable limits. E.g., using m* = 0.04mo 
instead of m* = 0.02m leads to (5(11) = 0.82 meV, 
(5(15) = 0.59 meV. 




FIG. S7. Structure of the bubble diagram corresponding to 
the analytic expression of Eq. (Sll). 
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